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What is the probabihty that a macroscopic void will spontaneously arise, at a specified time T, in an 
initially homogeneous gas? We address this question for diffusive lattice gases, and also determine 
the most probable density history leading to the void formation. We employ the macroscopic 
fluctuation theory by Bertini et al. and consider both annealed and quenched averaging procedures 
(the initial condition is allowed to fluctuate in the annealed setting). We show that in the annealed 
case the void formation probability is given by the equilibrium Boltzmann-Gibbs formula, so the 
probability is independent of T (and also of the void shape, as only the volume matters). In the 
quenched case, which is intrinsically non-equilibrium, we evaluate the void formation probability 
analytically for non-interacting random walkers and probe it numerically for the simple symmetric 
exclusion process. For voids that are small compared with the diffusion length VT, the equilibrium 
result for the void formation probability is recovered. We also re-derive our main results for non- 
interacting random walkers from an exact microscopic analysis. 
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I. INTRODUCTION 

In thermal equilibrium, fluctuations of macroscopic quantities are completely described in terms of free energy of the 
system f]\. Non-equilibrium fluctuations are much harder to study, especially when they are large. Anomalously large 
fluctuations, also called large deviations, of macroscopic observables in non- equilibrium steady states have recently 
attracted a lot of attention. Several lattice gas models have been investigated and it has been found that the 
distribution of fluctuations in non-equilibrium steady states, as described by the large deviation functional [2], can 
exhibit qualitatively new features, such as non- locality and phase transitions (see reviews [3, 4] and references therein). 

In this work we examine large deviations which are constrained to occur at a specified moment of time. We are 
aware of two one-dimensional settings where questions of this type have been already addressed. One of them deals 
with large fluctuations of integrated current for a step- like initial density profile p-Sl- The other deals with large 
deviations of dynamical activity, defined as the number of particle moves that have taken place over a given time 
window in finite systems with exclusion interaction [9\ . In general, problems of these type demand an explicit account 
of the system dynamics. 

Here we consider a different type of dynamic problem which allows us to study large deviations, for a whole class 
of lattice gases, in arbitrary spatial dimension. Starting with an infinite d-dimensional lattice filled with particles at 
constant density n, we study the probability of formation of a void, viz. a large empty region of a given shape with 
characteristic linear size L, at a specified time T. A remarkably similar characterization of large deviations has been 
considered in the context of quantum many-body systems, see e.g. 10, llj and references therein, under the name of 
"emptiness formation" . 

We will study the probability of void formation in two different settings — starting from a fluctuating equilibrium 
state (the so-called annealed setting) or from a deterministic state (the quenched setting). In the annealed setting 
the gas remains in equilibrium, for any T, in the process of void formation, as we explicitly show in this work. The 
quenched setting is intrinsically non-equilibrium, since it describes the void formation process as the gas is evolving 
on its way to equilibrium. 

Our analysis is based on a coarse-grained approach which directly probes the long-time limit and is applicable to a 
whole family of diffusive lattice gases of interacting particles. This approach is rooted in the macroscopic fluctuation 
theory (MFT) of Bertini, De Sole, Gabrielli, Jona-Lasinio, and Landim see also Refs. [1, H, The MFT is a 
generalization of the low-noise Freidlin-Wentzell theory [ij] which in turn is an analog of the WKB approximation 
of quantum mechanics. Similar approaches have been also developed for lattice gases where, in addition to diffusive 
transport, there are on-site reactions among particles p^l - [T7j . In its standard form, the MFT [l^l holds for lattice gases 
where the transport is unbiased and diffusion- like. In addition to non-interacting random walks (RWs), well-known 
examples of such lattice gases are the simple symmetric exclusion process (SSEP) 18-[25j. the Kipnis-Marchioro- 
Presutti (KMP) model |26l-[28|. and the symmetric zero range process (ZRP) [18|, |29|, |30[ . In the SSEP a particle can 
hop to a neighboring site if that site is empty; if it is occupied by another particle, the move is disallowed. The KMP 
model describes a lattice of mechanically uncoupled harmonic oscillators which randomly redistribute energy among 
neighbors. The ZRP describes interacting random walks: A particle at site i can hop to a neighboring site with rate 
a that depends on the number of particles ni on the departure site i. For these and other similar lattice gases, a 
hydrodynamic, or mean-field, description is provided by a diffusion equation 

dtp^W- [D{p)Wp] (1) 

for the average density p{r,t) [3l|. The diffusion coefficient D{p) is constant in the simplest models {e.g., for the 
RWs, the SSEP, and KMP), but generally it depends on the density; hence the diffusion equation JT]) is generally 
non-linear. 

At the level of MFT, fluctuating diffusive gases are fully characterized by D{p) and an additional function a{p) 
which describes equilibrium fluctuations [13) US- Table I lists the functions D{r) and a(r) for the four aforementioned 
models: the RWs, the SSEP, the KMP and the ZRP (in the latter case we assume that a{r) is a monotonically 
increasing function [sll). It also gives, for these models, the equilibrium free energy density F{r). The free energy 
density is related to D(r) and cr(r): 

cr(r) 

(Hereinafter, the prime denotes the derivative.) Equation ^ follows from the fluctuation-dissipation theorem [3l.[l8|. 
and it also naturally emerges from the MFT formalism j8j. 
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Model 


D{r) 


a(r) 


F(r) 


RWs 


1 


2r 


rlnr — r 


SSEP 


1 


2r(l - r) 


r In r + (1 — r) ln(l — r) 


KMP 


1 


4r-2 


-(1/2) Inr 


ZRP 


a'{r) 


2a(r) 


J' du In Q(ti) 



TABLE I: Functions D(r),a{r) and F{r) for non-interacting random walkers (RWs), the SSEP, the KMP and the ZRP models. 

In this paper we derive and analyze the MFT equations and boundary conditions describing formation of a void 
in a whole class of diffusive gases in arbitrary spatial dimension. We show that, for these systems, the probability 
V of void formation has a universal scaling form: InV ~ —T'^/'^S{L/y/T,n). In the annealed setting, the large 
deviation function S is such that the resulting V is independent of T and described by the classical Boltzmann-Gibbs 
equilibrium formula. Correspondingly, the optimal (most probable) time history of the gas density field in the process 
of void formation coincides in this case with a time-reversed history of the mean-field relaxation process, that is with 
a timc-rcversed solution of the diffusion equation. 

Finding S analytically in the quenched setting is a hard problem, since the optimal time history of the gas density in 
the process of void formation is different from the time-reversed solution of the diffusion equation. We have only been 
able to solve this problem analytically for non-interacting RWs. In general, S depends on the void shape. For short 
times (equivalently, large voids), the large deviation function is a genuinely non-equilibrium quantity that admits an 
integral representation with an interesting geometric flavor: 5 ~ n /y dX [I?(X)]^. Here V is obtained by rescaling 

all the coordinates of the void by the characteristic diffusion length v^iT, and I'(X) is the distance between the point 
X inside the rescaled void V and its boundary dV. In particular, this implies that the spherical void is the least 
probable among voids of the same volume. For long times (equivalently, small voids), the large deviation function S 
becomes shape-independent, and the equilibrium result for V is recovered. The equilibrium result is also obtained 
when maxPfX) <C 1. 

X 

Even in the relatively simple case of non-interacting RWs, the MFT formalism turns out to be quite instructive. 
Not only it accurately predicts the logarithm of the void formation probability, but it also gives the optimal time 
history of the gas density field. This history strongly depends on the ratio L/\/T. For interacting lattice gases, the 
void formation probability, and the optimal density history, can be found by solving the MFT equations numerically, 
as we demonstrate for the SSEP. 

The remainder of the paper is structured as follows. In section|TT]wc consider a diffusive gas with arbitrary D{q) and 
cr{q) and present the MFT equations and boundary conditions for the void formation problem; details of the derivations 
are given in Appendix A. In section IIIII we solve the void formation problem in the annealed setting. Sections IIVI 
and IVl deal with the quenched setting: Analytical results for the RWs are established in section HVl while numerical 
results for the SSEP are given in section |Vl Our main findings are briefiy discussed in section IVll In Appendix B we 
outline an exact microscopic theory of the void formation for the RWs, both in the quenched and annealed settings. 
This microscopic theory yields the void formation probability and the expected system configuration dX t — T which 
coincide, in the long time limit, with the corresponding results obtained from the MFT formalism. 

II. MACROSCOPIC FLUCTUATION THEORY OF VOID FORMATION 

A. Governing equations and boundary conditions 

Our analysis employs the MFT [T^ [l3j : a coarse-grained formalism which directly probes the long-time limit 
and is valid for a whole family of diffusive lattice gases. The MFT can be formulated as a classical Hamiltonian 
field theory, where the number density q{yi, t) plays the role of "coordinate" , and the conjugate field p(x, t) (which 
can be viewed as the magnitude of fluctuations) is the "momentum" . In Appendix A we present a derivation of the 
Hamilton equations and boundary conditions for the void formation problem in the annealed and quenched settings. 
The derivation starts from fluctuating hydrodynamics: a Langevin-type partial differential equation for q{'x.,t) whose 
deterministic part coincides with Eq. (P), whereas the (multiphcative) noise term includes a{q) [ll]. Being interested 
in the long-time behavior of the probability distribution, we arrive at a variational problem (see Appendix A) which 
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leads to two coupled partial differential equations for g(x, t) and p(x, i): 

dtq = V- [D{q)Vq - a{q)Vp] , (3) 
dtp = -D{q)V^p-]^a'{q){Vpf. (4) 

Equations ([3]) and (j4]) are Hamiltonian, since they can be written as 

dtq = SH/5p , dtp = -5H/Sq . (5) 

Here 

H[q{^,t),pi^,t)] = J d^H (6) 

is the Hamiltonian, 

H{q,p) = ~Diq)Vq ■ Vp + ^C7(g)(Vp)' . (7) 

The spatial integration in Eq. (jS]), as well as in a number of equations below, is over the entire space. For a given 
lattice gas model, specified by D{q) and cr(g), the same Eqs. dSj and (|4]) arise when one studies large deviations of 
different quantities in different settings. Boundary conditions in space and time, that complement Eqs. ([3]) and (j4|. 
vary from problem to problem. Here we study the formation of a void in an infinite system. Hence we demand that, 
at a specified time t — T, a. (simply connected) void of a given shape is observed: 

q{x, T) = inside the void. (8) 

We emphasize that we do not specify the density profile g(x, T) outside the void — it will emerge from the solution 
of the problem as the density profile that maximizes the probability to observe the void ([5]). 

The void formation probability not only depends on the dynamics of the underlying microscopic model during the 
time interval < t < T, but also on the initial condition. At the macroscopic level we want the initial density to be 
uniform. One way to achieve it is to start, in the microscopic formulation, with a deterministic constant density. The 
void formation probability is then obtained by averaging only over stochastic realizations of the dynamics over the 
time interval < t < T. Alternatively, we can allow equilibrium fluctuations of the initial condition and average over 
both these fluctuations and stochastic realizations of the dynamics. These two types of averaging are called quenched 
and annealed, respectively [11 . This terminology suggests an analogy with the quenched and annealed averaging in 
disordered systems, although there is no disorder in the present situation. Similarly to disordered systems, the analysis 
tends to be simpler in the annealed case, as we will see shortly. 

We now summarize the rest of boundary conditions for the MET equations ^ and As shown in Appendix A, 
the maximization of the void formation probability yields the following boundary condition at t — T: 

p(x, T) = outside the void, (9) 

for both annealed and quenched settings. Essentially, Eq. ^ states that fluctuations outside the void at i = T must 
vanish. The boundary condition at t = does depend on the setting. For the quenched setting the initial condition is 

q{x., 0) = n everywhere. (10) 

The annealed setting assumes equilibrium fluctuations in the initial condition, i.e., the density profile at t = is 
chosen from the equilibrium probability distribution corresponding to density n. For the annealed setting the initial 
condition for Eqs. ([3]) and (jlj, 

p(x,0)= J-'[g(x,0)], (11) 

establishes a relation between the most probable initial density profile and the corresponding profile of p(x, 0), see 
Appendix A. The function J-{r) is simply related to the free energy density F{r), differing from it by a linear function. 
Namely, T{r) obeys Eq. ^ and the additional relation T'{n) — 0. For concreteness, we also set F{n) — 0, so that 

^(,.).r.,/%c^. (12, 

Jn Jn Cr(C) 



5 



Finally, the boundary conditions at x — ^ oo are 

q{x oo,t) — n, p(x ^ cx), = 0. (13) 

Although Eqs. ([3]) and ^ and the boundary conditions represent a complete set for the quenched and annealed 
settings, there is an important corollary [valid under certain conditions on the functions D(r) and (j{r)] that can be 
very useful for solving the problem. This corollary is 

p(x, T) = — oo inside the void, (14) 

in both quenched and annealed settings, see Appendix A. 

One additional comment is in order about the fields q and p at t = in the quenched setting. Although the initial 
density profile (/(x, 0) = n is fiat here, the (a priori unknown) initial momentum p(x, 0) is nonzero. Therefore, the 
initial state, as described by the MFT, is non-deterministic. How to reconcile this fact with our definition of the 
quenched setting as the one starting from a deterministic density profile? The solution comes from the realization 
that, for the MFT formalism to become valid, one should wait for a time which is much longer than the characteristic 
microscopic time of the system: the time of a single particle move. During this time (which is assumed to be very 
short compared to the macroscopic time scales that the MFT formalism can only deal with) the deterministic density 
profile rapidly evolves into a fluctuating profile with constant average density n. In contrast to the annealed setting, 
however, these fluctuations are not in equilibrium: they are determined by the a priori unknown held p(x, 0) that 
depends on £ and is different from the equilibrium p-field described by Eq. (jlip . 

Solutions of Eqs. ^ and (Bl) with vanishing momentum, g(x, t) — p{x., t) and p(x, t) — 0, are deterministic; they are 
called relaxation solutions [3J| . For the relaxation solutions Eq. @ is satisfied identically, whereas Eq. 1^ reduces to 
the diffusion equation ([T|). Solutions with p{x, t) ^ are called activation solutions, for these solutions q{x, t) ^ p{x, t). 
The void formation obviously demands an activation solution. Once the activation solution g(x, t) and p(x, t), obeying 
the boundary conditions, is found, we can evaluate the action of the Hamiltonian system and (U): 

S^JdxJ^ dt {pdtq - n) = ^ J dt<j{q){\/pf. (15) 

For the quenched setting, this action yields the void formation probability: InT' ~ —5. For the annealed setting, one 
also has to account for the "cost" 5*0 of creating the optimal initial density profile (/(x, 0) in the equilibrium gas of 
average density n. This cost is given by the Boltzmann-Gibbs equilibrium formula, so 

Sa = j dx^[g(x,0)] (16) 

(recall that we defined T{r) so that J^{n) = 0.) Therefore, for the annealed setting we have InT' ~ —(5*0 + S). 



B. Dynamic scaling of the void formation probability 

For a given diffusive lattice gas, the action S depends on the characteristic size of the void i, the formation time 
T, and the gas density n. Let us rescale time by T, t/T — > t, and the distances by the diffusion length, x/-\/T x. 
Equations (jS]) and Q remain invariant under this transformation. The boundary conditions ([5]) and ^ [and the 
corollary (jl4p ] remain the same except that T is replaced by 1, and the physical void is replaced by the void of rescaled 
size L/^/T. The rest of the boundary conditions do not change. Equation ((T5|) becomes 

S{L,T,n)^^T''^' J dx j\ta{q){Vp)\ (17) 

and Eq. (1161) also acquires factor T^l"^ . Now, parameters L and T only enter the rescaled problem via the combination 
L/\/T. Therefore p{x,t) and q{x,t) may depend on L and T only via this combination. This implies a scaling form 

S{L,T,n) = T'^/^S{L/VT,n) (18) 

in both quenched and annealed settings, and for a whole class of symmetric diffusive lattice gases. The problem, 
therefore, is reduced to finding the large deviation function S{L/y/T, n) which coincides with the mechanical action of 
the rescaled problem. In the next section we show that in the annealed setting the large deviation function drastically 
simplifies, viz. the sum S'o + iS' is given by the Boltzmann-Gibbs equilibrium formula, so the void formation probability 
is independent of T and of the void shape. 
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Model 


In ^annealed 


RWs 




SSEP 


-V^ln-r^ 

1 — n 


KMP 


— V X oo = — oo 


ZRP 


— V [nlnn — Jj" du In a{u)] 



TABLE II: InV in the annealed setting for the RWs, the SSEP, the KMP, and the ZRP. The results depend only on the average 
gas density n and the volume of the void V. 



III. VOID FORMATION IN ANNEALED SETTING 



A. General 



The annealed setting turns out to be simple. The reason is that the initial condition belongs to the invariant 
equilibrium manifold of Eqs. ^ and described by the local relation 

p(x,t) - J-'[g(x,t)] (19) 

between q and p P. Using this relation along with Eq. ([2]), one can see by direct calculation that q{x,t) is governed 
by the time-reversed deterministic equation 

dtq = -'7 ■[D{q)SIql (20) 

as expected for reversible models, like ours, in equilibrium. Using Eqs. ([U, ([9]) and (fTO]). we obtain the full density 
profile dX t = T: 

g(x,r) = ^' "^«ide the void, ^^1) 
n, outside the void. 



With this condition we can solve the "anti-diffusion" equation \2Q\ backward in time. In this way we obtain the 
optimal density history of the gas, whereas q(x,0) yields the optimal initial condition. Now we can calculate p{x,t) 
from the local relation and determine S and Sq from Eqs. and ((TC)) . The following shortcut, however, makes 
these calculations redundant. The creation of optimal initial profile g(x, 0) at t = 0, followed by the void formation at 
t = T, can be described as a single extended activation trajectory q{'x.,t) that starts, at t = — oo, from the flat state 
q{x^t = — oo) = n, acquires the optimal shape (/(x, 0) at i = and ends by forming the desired void at t = T. This 
extended trajectory belongs to the invariant equilibrium manifold. Therefore, the cost of creating a void is determined 
by the Boltzmann-Gibbs equilibrium formula, and we obtain 

InT'anncalcd " d^J^[q{x,T)] ^ - J dxJ-(O) = -J"(0)F, (22) 

void 

where V is the volume of the void, and we have used our convention -Fi^n) = 0, see Eq. (jl2p . The equilibrium result 
(l22t is independent of T and of the void shape. 

Table II yields In 7^ for the models listed in Table HI Note that In 7^ is proportional to n for the RWs reflecting their 
non-interacting character. For the KMP model, J^(0) = oo which implies zero void formation probability; this is also 
evident from the definition of the microscopic model (26)]. For the SSEP, J^(0) diverges as n — J' 1, again as expected 
from the microscopic model. 

Although the probability V depends, via J-{0), on both D{r) and a-{r), the density histories only depend on D{r) 
but not on cr{r). In other words, all diffusive lattice gases with the same D{r) and the same average density have 
identical optimal density histories of void formation in the annealed setting. We now present more details on the void 
formation in one dimension. 



B. Constant diffusion coefficient 



Let D{q) = 1. Such a density-independent diffusion coefficient characterizes, e.g., the RWs, the SSEP, and the 
KMP. In one dimension, the void is a segment (we always tacitly assume that the void is a connected set). We can 
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x=x/V4r 

FIG. 1: (Color online) 

Void formation in the annealed case for models with D = 1. Shown are the optimal density profile histories q{X,t)/n at times 
t/T = 0, 0.9, 0.99 and 1 for ^ = L/^/W = 0.05 (left panel), at times t/T = 0, 0.5, 0.75 and 1 for ^ = 0.5 (right panel), and at 
times t/T = 0, 0.75 and 1 for £ — 5 (bottom panel). The activation solutions are time-reversed relaxation solutions. 



set the void to be the [~L, L] segment. Thus Eq. ([2T|) becomes 

q{x,T)=nH{\X\-e)=nH{\x\-L), (23) 

where H(z) is the Heaviside step function. The optimal density history of void formation is described by the solution 
of the hnear anti-diffusion equation: 

r .^ n f £-X \ n f i + X \ 
q{x, i) = - erfc , + - erfc , , (24 

2 \^T^J 2 \^T^J' 

where X = a;/\/4T, £ = L/a/4T and erfc(z) ~ (2/y^) exp(— is the complementary error function. The 
optimal initial density profile, therefore, is 

72 Tl 

q{x, 0) = - erfc {I ~ X) + - erfc {I + X) . (25) 

Some examples of the density history are shown on Fig. [1] The left panel corresponds to small I. In this situation, 
the optimal initial profile q{x, 0) is almost flat. As a result, most of the actual void formation occurs towards the 
end of the time interval Q < t < T . The bottom panel shows the opposite regime (large t) . Here the optimal initial 
density proflle already has a pronounced dip: the equilibrium fluctuations had to do most of the job already at t < 0. 
Equations ((22)) and ([25)1 for the RWs also follow from exact results of the microscopic theory (Appendix B). 



C. Nonlinear diffusion 



For a density-dependent diffusion coefficient, the optimal density histories [described by the nonlinear anti-diffusion 
equation (j20l) with the boundary condition (j23l) ] are more interesting. One fascinating phenomenon appears in a class 
of models satisfying D{Q) = 0. In such models, in the process of void formation, the gas density vanishes in a finite 
region of space already at an earlier time t < T. The empty region grows with time until the complete void ([25]) is 
formed at t = T. For sufficiently small Ljs/T one has r > 0, and the initial density q(x, 0) is everywhere positive. 
Otherwise the empty region is already present in the optimal initial density (/(x, 0). 
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We now give more details in the case when the diffusion coefScient has a simple algebraic form: D{q) = with 
fc > 0. When L/^/T < ^fc, the initial density is everywhere positive. If L/^/T > ^k, the initial density q{x, 0) vanishes 
on a finite interval. The factor comes from the self-similar solution, 

q,{x,t)^n^(^-^^), (26) 

of an auxiliary relaxation problem: 

dtq = d.,{q''d.,q) (27) 

on the interval < oo, with the initial condition being a step function: q{x, 0) = nH(—x). A distinctive feature of this 
class of problems (see e.g. [SA-Sl]) is the (semi-)compact support: there exists a finite point in space, Xk{t) = £,k Vt, 
so that the solution qs{x, t) is positive at a; < Xk{t) and zero at x > Xk{t). The shape function 0(^) solves the ordinary 
differential equation 



(1/2) ^(^' = (28) 



with boundary conditions 0(^ — —oo) — 1 and 0(^ — )■ oo) = 0. This problem can be easily solved numerically by a 
shooting method. Figure [5] gives an example of numerical solution for fc = 1, that is D{q) — q. Here ^k=i — 1-239. 

Now consider Eq. (|27l) when the initial condition includes two step functions, as described by the right hand side 
of Eq. ([23| . Because of the semi-compact support of the similarity solution (|26|). the solution of this problem, at 
sufficiently short times t < r, is a sum of two counter-propagating similarity solutions of the type ([26]): 



+"'^(;^^2^) ' *<^- (29) 

The character of solution changes when the edge points Xj7 = —L + £,kVt and x'^ = L ~ £,kVt of the two similarity 
solutions meet. This occurs a.t t = L'^/£,k- As the activation trajectory is a time-reversed relaxation trajectory, one 
can easily obtain the condition L/^/T < for the positiveness of q{x, 0) everywhere, and L/^/T > for the presence 
of an empty region already at i = 0. Again, the same density histories will be observed for all models with the same 
D{q), independently of cr(g). 

Figure [3] shows a complete density history for L = 1 and T = 2, obtained by solving numerically the anti-diffusion 
equation (|20|) . for fc = 1, with the boundary condition ((23)) . As L/^/T < in this example, the initial density is 
everywhere positive. An empty interval appears at time t — T — L^^k=i — 1-"^"^ ^ corner in the density profile. 
Finally, InPanncaicd — — 2L X 2n = —2L x 2 x 1 = — 4_L in this case, see Table 2. For comparison, InPanncaiod — —2L 
for the RWs with the same average density n = 1. 



IV. VOID FORMATION IN QUENCHED SETTING: NON-INTERACTING RANDOM WALKERS 



A. General 



In the quenched setting one needs to solve Eqs. (|31) and (|1]) with the boundary conditions at t = and T described 
by Eqs. (jlOp and (fH|) . Now the system does not lie in the invariant equilibrium manifold ([T^ . and the problem 
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X 

FIG. 3: (Color online) 

Void formation in the annealed case for models with D{q) — q. The parameters are n = 1, L = 1 and T — 2. Shown is the 
optimal density history q{x,t) at times t = 0, 1, 1.35, 1.9 and 2. At time t ~ 1.35 an empty interval appears. The activation 

solutions are time-reversed relaxation solutions. 



of void formation appears intractable for a generic lattice gas. We have succeeded in solving it analytically only for 
RWs, where Eqs. ([3]) and become 



dtq - V-(Vg-2gVp), 
dtp = -V2p-(Vp)2. 



(30) 
(31) 



The solution can be obtained via the Hopf-Cole canonical transformation defined by the relations Q = qe ^ and 
P — e^, see e.g. [l^. The generating function of this transformation can be chosen as / dx^{q,P) = J dxqlnP. 
The new Hamiltonian is J dxH, with density H = —VQ ■ VP. The Hamilton equations become 



dtQ = V^Q, 

dtP -- 



V^P. 



(32) 
(33) 



The mechanical action along an activation trajectory can be written as 



S = J dyi dt (pdtq-H) 

= J d^J^ dt (PdtQ-'fdj + j dx^{q,P) 
^ dt J dyi{P\/^Q + \/Q -S/P) + J dx^q,P) 



(34) 



The integral over x in the first term vanishes by virtue of Green's first identity and the boundary conditions at x — >■ oo. 
As a result, the action 

S ^ J dii<i>{q,P) = y"dx [g(x,r) lnP(x,r) -g(x,0) lnP(x,0)] (35) 
is fully determined by the initial {t = 0) and final {t — T) states. In view of Eqs. (j9]) and (fT4|) we have 



inside the void, 

1 outside the void. 



p(x,r) = 

Using Eqs. ([5]) and we can reduce Eq. ([55t to 

S = -n I dx lnP(x,0) 



(36) 



(37) 



Now, P(x, 0) can be easily found by solving the linear anti-diffusion equation (|33l) with the boundary condition (|36p . 
Let us first consider the one-dimensional case. 
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B. RWs in one dimension 



In one dimension, the boundary condition 
with this boundary condition, we find 



36| becomes P{x,T) = H{\x\ — L). Solving the anti-diffusion equation 



At t = this yields 



P(x, 0) = ^ erfc - X) + i crfc {(. + X) . 



Plugging it into Eq. (|37| . we arrive at 



\uV{L,T,n)~ S = nV4Ts(£) 



where 



/CO 
dX In 
-OO 



erfc(£ -X) + erfc(£ + X) 



(38) 



(39) 



(40) 



(41) 



One can see that h\V{L,T,n) indeed exhibits dynamic scaling, as expected for all diffusive gases, see Sec. Ill Bl The 
linear dependence on the density, S oc n, reflects the non-interacting character of RWs. 
Let us find asymptotic behaviors oi s{£). For £ Eq. (HIT) yields 



&{£) = 2£+\/-> 



so that 



lnP(L,T,n) 



n 2L 



(42) 



(43) 



We emphasize two important features. First, the leading-order term, lnP(L,r) ^ —2nL, coincides with the annealed 
result for the RWs, see Table 2 with V — 2L. This is anticipated as in the limit of £ — > (small L or large T), the 
system is given sufficient time to approach equilibrium and exploit equilibrium fiuctuation for creating an optimal 
initial density profile which facilitates the void formation. We expect this feature to hold for interacting diffusive 
gases and confirm this expectation numerically for the SSEP in Sec. |Vl Second, the void formation probability in the 
quenched setting is (exponentially) smaller than the void formation probability in the annealed setting. This feature 
is general: by fixing g(x, 0) = n, we narrow the class of density variations in the problem of maximizing V. As a 
result, Pquenched < ^annealed, See also Appendix A2. 

In the far-from-equilibrium limit, £ ^ 1, we can employ the large- ?/ asymptotic of erfc(?/). 



rfc(y) = 



1 



(44) 



to extract the asymptotic of the integral in (|41l) . The main contribution to the integral is gathered in the region 
\Z\ < £, and we obtain 



5{£) = -£3 + 2£\ii£- 



This leads to 



\nr{L,T, n) 



[ln{4Tr) ~2]£ + 0{l). 



6T 



2L In 



(45) 



(46) 



Note that the leading term of the asymptotic pS)) . (2/3) ^'^ comes from approximating the integral in Eq. (|^T|) as 

ill dz{z — £)^. This is the integral of the squared distance from the boundary of the rescaled void over the rescaled 
void. As we will see shortly, this geometric property holds, at ^ ^ 1, for RWs in any dimension. Figure 2] depicts 
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I 

FIG. 4: (Color online) 

Function s{£) from Eq. I|41|) (solid line), the small-^ asymptotic (|42|l (dashed line) and the large-^ asymptotic (|46p (dotted 

line) . 



5{i) along with its asymptotic behaviors. One can see that it is exponentially less probable that a void of a given size 
appears in a short time, than in a long time. 

Now let us see how the optimal density field evolves in time. Using Eq. ([39]) and the condition q{x, 0) = n, we 
obtain 



Qix,0) 



2n 



Pix,0) erk {£ - X) + evfc {£ + X) 
Solving the diffusion equation p2p with this boundary condition yields 



AT 

Q(x,0< t < T) = n\ — 
V Tri 



; erfc (^ — z) + erfc + z) 
Equations (|38p and (|48| together with relation q = QP give us the optimal density. In particular, 

2nH{\x\-L) r e-(^-^)'d2 



q{x,T) 



erfc {t- z) + erfc (£ + z) ' 



(47) 



(48) 



(49) 



The optimal density history in the quenched case does not coincide with the time-reversed solution of the diffusion 
equation. This can be seen, for different values of from Fig. [S] For ^ > 1 there is a striking difference between 
the annealed and quenched density histories outside the void. For / ^ 1 (the left panel), the annealed and quenched 
optimal density histories become similar. In particular, most of the void formation in this case occurs towards the 
end of the time interval < t < T, when the system can exploit (almost) equilibrium fluctuations. 

As one can see, the gas density profile outside of the quenched void at t = T is quite non-trivial. For instance, there 
are density cusps at the void boundaries x ~ ±L. The magnitude of the density at the peaks is 



(?(±L,T) = ^ / dz 

Vl" J-oo 



erfc(£ — z) + erfc(£ -|- z) 



(50) 



For £ < 1 Eq. (gH]) becomes 



q{\x\ >L,T)~n{ 1 + J-ee-^ 



which is close to n everywhere. On the contrary, for £ 3> 1 the gas density at the boundaries is much greater than n: 

q{±L, T) ~ n _^ In £ -f i ln(47r) + 0{e-^)] . (51) 

It rapidly decays, however, as a function of x, and approaches the initial density n. The decay occurs in a boundary 
layer of width ^ 1. Indeed, writing \X\ — £ — £^/£ and taking the limit £ — > oo with kept finite, we can simplify 
the integral in Eq. (|49p and obtain 



g(x,T) = n£2^(0, i^iO 



1 



-2« _ 



-2« 



(52) 
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FIG. 5: (Color online) 

Void formation in the quenched case for independent random walkers. Shown are the optimal density histories q(X, t) at 
times t/T = 0, 0.9, 0.99 and 1 for ^ = 0.05 (left panel), at times t/T = 0, 0.5, 0.75 and 1 for ^ = 0.5 (right panel), and at times 

t/T = 0, 1/3, 2/3 and 1 for ^ = 3 (bottom panel). 



In the region of 1 <C ^ ^ or equivalently ^ ^ <C \X\ — ^ ^ 1, the gas density exhibits a power-law decay in x: 

n 2nT 



(53) 



2{\x\-er {\x\^L)^- 

Equations (PU)) . (PT|) and (|49l) can be also obtained from exact resuhs of the microscopic theory, see Appendix B. 



C. RWs in higher dimensions 



Now let us return to Eq. ([37]) and consider formation of a void of any simply-connected shape in d dimensions. To 
calculate In 7^, we need to solve the anti-diffusion equation (1331) backward in time with the boundary condition (j36l) . 
evaluate the result at t = and plug it into Eq. (|37p . In this way we obtain 



In-P 



-n(4T)'^/2srf, 



where 



Sd = — / rfZ In 



1 - TT-"^^ / rfYe-(2-Y) 
Jv 



(54) 



(55) 



The integration over Z is performed over the whole c?-dimensional space, whereas the integration over Y is performed 
over region V which is obtained by rescaling all coordinates of the physical void by •\/4T. 



1. Small voids 



Consider voids with all characteristic dimensions much less than the diffusion length v 4T. Such voids have a small 
rescaled volume, |V| ^ 1. Expanding the logarithm in (j55p and keeping the leading term we get 



(56) 
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Exchanging the order of integration we find that the integral reduces to |V|. This leads to 

Vd{V,T) 0^ e-""" , (57) 

which coincides with the annealed void formation probability for RWs, see Table II. Expanding the logarithm in ([55]) 
to higher orders, we find 

Srf(V) = |V| + i(2^)-'^/2|Vp + ... (58) 

Surprisingly, the sub- leading term also depends only on the rescaled volume |V| and is independent of the shape of 
the void. In one dimension |V| — 2£, and Eq. (|58|) coincides with Eq. (j42|) . 

The equilibrium result (j57p actually holds under a much weaker assumption on the spatial dimensions of the rescaled 
void. It suffices to demand the following strong inequality: 

maxX>(X) < 1, (59) 

where I'(X) is the shortest distance from X to the boundary dV of the rescaled void. Under condition ([59)) . the field 
P(X, 0) that determines the action ([37]) is close to unity in all points inside the void, validating the approximation 
([5B)) leading to Eq. ([57]). 



2. Large voids 

Now consider voids with all characteristic dimensions much greater than the diffusion length 
a large rescaled volume, |V| ^ 1. In this case it is convenient to rewrite ([55]) as 



Srf(V) 



dZ In 



-d/2 



dYe 



-(Z-Y) 



4T. Such voids have 



(60) 



where the second integration now goes over the exterior V = M'' — V of the rescaled void V. The term inside the 
logarithm in Eq. ()60p is very close to 1 when Z S V, apart from the region very close to the boundary of the rescaled 
void V, but the contribution from this region is very small as the rescaled void is large. For the points Z e V, the 
term inside the logarithm is very small. This provides the dominant contribution to the first integral in (j60p and 
implies a simplification. Indeed, let Ym be the point on the boundary dV of the rescaled void V which is closest to 
Z. In other words, I'(Z) = |Z — Yz| is the shortest distance from Z to the rescaled void boundary dV. Then a crude 
estimate of the internal integral in Eq. (|60p , 



Jv 



suffices to provide a correct leading contribution to Eq. ([50]) . Thus we arrive at 



Sd(V) 



dZ [V{Z)f 



(61) 



(62) 



a purely geometric result. 

Now we are in a position to ask the following question. Given that a void of volume V has formed in a diffusive 
lattice gas at time T, what is the most likely void shape? To get insight, let us first apply Eq. ([62]) to a d-dimensional 
spherical void of rescaled radius I = R/^/^T 3> 1. We obtain 



,ball 



dvd I dr r' 
Jo 



4-1 



'd+2 



2«7'/''|V|H 



(d+l)(d + 2) (d+l)(d+2)' 



(63) 



where — 7r'^/^/r(l + d/2) is the volume of the unit sphere in d dimensions. Combining Eqs. ([5^ and ([M| and 
recalling that |V| = Vdi"^ we obtain 



Vd{V,T,n)^exp 



V 



(d + l)(d + 2) 2T 



(64) 
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Similarly, for a cube-shaped void we find 

cube 11^1^^^ 

2(d+l)(d + 2)- ^^^^ 

An inspection shows that, for a given |V|, and for d > 1, 5'^^^° < sj^**". That is, at short times, the formation of a 
spherical void is (exponentially) less probable than the formation of a cube-shaped void of the same volume. But the 
cube is also far from the most probable shape at a given void volume. This becomes clear when one uses Eq. (|62p 
to calculate Sd for a d-dimensional rectangular parallelepiped (cuboid). For example, for a two-dimensional rectangle 
with rescaled sides of lengths a ^ 1 and 6 ^ 1, a > 6, we obtain 



rectangle _ b^C^a - b) _ \V\^ 2k - 1 

^2 24 " '24 ^^^> 



where k — a/b > 1 is the aspect ratio of the rectangle, and |V| = ab is the rectangle area. The function k^^{2k ~ 1), 
entering Eq. (|66p. has its maximum at k = 1 (that is, for a square-shaped void) and is monotone decreasing at k > 1. 
Hence the minimum of Sd at fixed |V| is achieved at very large n (formally at k oo), when the rectangle becomes 
a long and thin filament. For thin filaments, however, the large- void approximation, Eq. (|62p. breaks down. On the 
other hand, in the limit of a 3> 1 and 6^1 the equilibrium asymptotic (|58p holds, which is independent, in the leading 
and sub-leading order, of the void shape. Generalizing to d dimensions, we see that the most likely void shape, at a 
given volume, is such that at least one characteristic dimension is much smaller than the diffusion length \/AT. In this 
case the void formation probability, in the leading order, is close to the equilibrium probability which is independent 
of the shape. At this level of accuracy we cannot distinguish between a whole variety of shapes: for example, between 
a pancake and a filament of the same volume in Sd. To determine the most likely shape, one may need to account for 
non-equilibrium sub-leading terms in that can still be important, as they appear in the exponent. 



3. Spherical voids 

Among voids of a given large volume, spherical voids arise with the lowest probability; this assertion is rather 
evident, although we do not have a rigorous proof. Let V be the rescaled ball of rescaled radius £. Writing |Z| = r 
and using spherical symmetry we get dZ = dvdr'^~^dr. Simplifying the integral in (|55p . we obtain 



Sd{e)^-dvd drr^-' \nfd{r,£) 



(67) 



where 



/,(r,£) = l-^-^/2e-'-' 



dYe 



2Z-Y-|Y|- 



|Y|<£ 



In three dimensions, for instance. 



S3ii) = -4n drr'\nh{r,l). 



(68) 



and f'i{r,t} simplifies to 



erfc(^ — r) + erfc(^ -I- r) e 



2y/TTr 



(69) 



Now we can return to the £ ^ oo limit and calculate a sub-leading correction to the already known leading asymptotic 
Using Eq. dH]), we obtain 



.3m = |^'.f^= 



ln£+i ln(47r)-if 
2 b 



(70) 
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V. VOID FORMATION IN QUENCHED SETTING: SSEP 



Now we briefly consider the void formation in the SSEP. For the quenched setting, Eqs. ^ and (HJ cannot be 
solved analytically. We solved these equations numericaUy in one spatial dimension, by using an iteration algorithm 
originally developed by Chernykh and Stepanov [s^ for evaluating the probability distribution of large negative 
velocity gradients in the Burgers turbulence. Different modifications of this algorithm have been used for evaluating 
large deviation functions of several lattice gas models, with and without on-site reactions [1, [l^ [l^. This algorithm 
is ideally suitable for the void formation problem, as this problem involves mixed boundary conditions in time: one 
condition for q, q{x, 0) — n, and another for p, see Eqs. ^ and The algorithm iterates the diffusion-type Eq. ^ 
forward in time and the anti-diffusion- type Eq. Q backward in time (see fsl for details). To suppress numerical 
instability, we used a linear combination of values of q and p from two previous iterations when solving for p and q, 
respectively [1, [s^ . The only significant difference between our present implementation of the algorithm and that of 
Ref. Q comes from the boundary condition p(| a; | < L,T) = —oo which is inconvenient numerically. To overcome this 
inconvenience, we performed the Hopf-Cole canonical transformation from q and p to Q = qe~^ and P = e^. The 
transformed Hamiltonian density is 



and the new Hamilton equations read 



As a result, the action becomes 



H^-d.,Qd^P-Q\d^,Py, (71) 



dtQ = d^^Q + 2d^iQ^d,P), (72) 
dtP = -d,,P + 2Qid,Pf. (73) 



S = J J [pdtQ-Hj dxdt + J ^{q,P)dx 

^ - j j Q'^idxPf dxdt~n j dxh\P{x,Q), (74) 

where $(g, P) — glnP, see Sec. IIVI A. The boundary conditions in time become 0)P(a;, 0) — n and P{x,T) ~ 
-ff(|a;| — L) — 1. The boundary conditions in x do not change: Q(|a;| ^ oo,t) = 1 and P{\x\ ^ oo,t) — 0. 




FIG. 6: (Color online) Void formation in the quenched case for the SSEP. Shown are numerically computed optimal density 
histories q{x,t) for L — 0.25 and T — 1 (left panel) and for L = 3 and T = 1 (right panel) at time moments t = 0, 1/3, 2/3 
and 1. The initial density n = 0.4 



We implemented this algorithm in Mathematical working with finite systems of reasonable sizes. The step-function 
entering the boundary conditions for P at t = T was smoothed a bit. The iterations converge rapidly. Mass 
conservation in the numerical box was used to monitor the accuracy. Having computed Q{x,t) and P{x,t), we 
determined q{x,t) = Q{x,t)P{x,t) and computed the action S by numerically evaluating the integrals in Eq. (|74p . 

Figure [5] shows two examples of numerically found optimal density histories q{x,t) for n = 0.4 and T = 1: for a 
small void, L = 0.25, that \s t — 0.125 (the left panel) and a relatively large void, L = 3, that is ^ = 1.5. One can see 
that for L — 0.25 the density history is similar to those for RWs with small I, as shown in the two upper panels of Fig. 
[SJ For large I this similarity must break down, as the gas density in the SSEP cannot exceed g = 1. Furthermore, for 
very large voids, £ ^ 1, the characteristic decay length of the density outside the void must behave as £, the same as 
the displaced mass. (This should be contrasted with the RWs where, at n = 1, the maximum gas density is ^ 1, 



16 



CO 




FIG. 7: (Color online). S{£,n) versus £ for the quenched setting of the SSEP with n = 0.4. Circles: numerical results obtained 
with the iteration algorithm. Dashed line: the annealed result 5 = 4£ ln[l/(l — n)], see the text. 

and the density decay length is 1/f <C 1, see Sec. IIVI B.) These two features of the SSEP can be discerned already at 
£ = 1.5, see the right panel of Fig. [S) We could not probe £ > 1.5, as the numerical accuracy became insufficient. 

Figure [7] depicts the ^-dependence of the rescaled action S{£, n) that we found numerically in a moderate range 
of ^ = Ll\/^T for n ~ 0.4. As expected, at small £ the rescaled action approaches that for the annealed case, 
iSanneaicd ~ 4i?ln[l/(l — n)\. The latter relation follows from Eq. (P^ . see Table II, and the dynamic scaling relation 
S{L,T^n) — T^^'^S{£,n), see Sec. Ill Bl As also expected, the rescaled action for the quenched setting is greater than 
the one for the annealed setting, so that Pquenched < 'Panneaied- 



VI. DISCUSSION 



We investigated the probability of macroscopic void formation in a class of diffusive lattice gases whose hydro- 
dynamic description is provided by the diffusion equation ([T]). The formalism of macroscopic fluctuation theory is 
perfectly suitable for the analysis of the void formation problem — it predicts the dynamic scaling behavior (|18|) . and 
it also yields the most likely density history of the system in the process of void formation. In the annealed setting the 
void formation probability turns out to be independent of T and given by the Gibbs-Boltzmann equilibrium formula. 
The quenched setting is harder to study, and we were only able to solve analytically the case of non-interacting 
random walkers. Fortunately, a relatively straightforward numerical treatment is feasible for interacting particles, as 
we have demonstrated for the simple symmetric exclusion process. An interesting avenue for the future research is 
to develop an analytical theory in the limit of large voids. The hope is to use asymptotic methods to circumvent 
the challenge of (most likely, inaccessible) exact solution of the MFT equations for interacting lattice gases. More 
specifically, it would be interesting to find whether the super- Gaussian tail InT' ~ — f{n)L^ /T, in analogy with our 
RW result (|46p. is universal for a whole class of interacting lattice gases. A similar question about the super-Gaussian 
tail ln7'(J) ~ ~J^/T of the probability to observe a very large integrated current J in one dimension was raised in 
and the arguments in favor of this behavior were given in Refs. HQ. 

We argued that the probability of formation of macroscopic voids is a meaningful way to characterize large deviations 
in classical diffusive lattice gases. As we already mentioned in the Introduction, a similar characterization of large 
deviations has been employed for quantum many-body systems, see e.g. [lol. ITllj. where it is known under the name 
of "emptiness formation" . One approach to finding the emptiness formation probability relies on exact solutions of 
integrable models [S^, and subsequent computation of the large-distance asymptotic behaviors, see [ll|. Another 
app roach is via an effective nonlinear hydrodynamic description that directly probes the large-distance asymptotics 
[lol[Tl| . In the latter description the probability of void formation corresponds to an activation trajectory: an instanton 
solution of a classical hydrodynamics (or rather "anti-hydrodynamics" where, in Id, the usual pair of acoustic waves 
in a compressible gas give way to a pair of aperiodic modes: one growing, the other decaying). On a qualitative level, 
this "anti-hydrodynamics" may be compared to the "anti-diffusion" that appears in the lattice gas settings we dealt 
with here. This analogy opens new exciting directions for future work. 
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Appendix A: Derivation of the MFT equations for void formation 

Here we derive the MFT equations and boundary conditions for void formation in both annealed and quenched 
setting. Rather than immediately focusing on the void formation, it proves useful to start with a more general problem 
and to demand qipi, T) — k(x) in a bounded domain $7; the void problem is then obtained by identifying Q with the 
void and setting k(x) — 0. 



For the diffusive lattice gases, the Langevin equation is [18 1 



dtq^V-[Diq)Vq+yMQ)^ix,t)\. (Al) 
Here ^{x,t) is a zero-average Gaussian noise, which is delta-correlated in space and in time, 

(6(x,i)0(x',i')) =^..'5(x-x')5(t-0, (A2) 

and the brackets denote ensemble averaging. We assume that q{\x.\ — oo, 0) ~ n = const. It proves useful to introduce 
the particle displacement u defined via 

g(x,t) = n-K V • u(x,t) (A3) 

and satisfying the boundary condition 

u(|x| ^ cx),0) = 0. (A4) 
Combining (jAip and jASj, we deduce a Langevin equation for u: 



dtU-D{q)Wq=^/^^. (A5) 
The probability of observing a large deviation of the displacement u(x, t) is, up to a pre-exponent, 

VH.,t)] ^ exp|-£rf. / I^^^i^^M^I . (A6) 



Here and in the following we assume that q is related to u according to Eq. (|A3p . 

1. Annealed setting 

Let us first consider the annealed setting. At t = 0, one starts from an (a priori unknown) optimal initial density 
profile g(x, 0): a certain realization of equilibrium fluctuations of the gas with average density n = const. The 
probability to observe this realization is given by the Boltzmann-Gibbs distribution and can be expressed via the 
function J-{q) defined in Eq. (|12l) : 



7'[g(x,0)] ^expj-y rfx^[g(x,0)]| , 



The joint probability to observe a large deviation of the displacement u(x, t) when starting, at t = 0, from q{x, 0) is 

V[u{x,t);q{x,0)]^exp{~S), 

where 

J Mq) 



I dt I ^ , dxT[q{xM. (A7) 
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Our task is to minimize S under condition that q{x,T) = k(x) in a bounded domain C R''. Let us define a scalar 
field p(x, t) (which will play the role of momentum) as the unique solution of the following Poisson-type equation: 



with the boundary condition 



V • (aVp) = V • (D Vq) - dtq 



p{\x\ -> oo,i) 0. 



(AS) 



(A9) 



As q tends to n = const at x — > oo, the right- hand- side of Eq. (jA8l) (the "charge density" of the Poisson equation) 
vanishes at x — > oo, and so the solution for p is indeed unique. Equation (|A8p coincides with the first of the MET 
equations, Eq. ([3]). 

An immediate corollary from Eqs. (jASp and (jA3l) is the relation 

(AlO) 



dtU — D Wq = — (T Vp + uj . 



where a; is an arbitrary solenoidal vector field, W ■ u) — 0, which goes to zero sufficiently fast at x — oo. Although 
UJ does not affect q{x,t), it does affect the action S. Indeed, upon substituting Eq. (|A10p into Eq. (jA7[) . the integral 
over X in the first term can be written as 



dx 



a{Vpf 



2^ 



dx V • {pLj) 



(All) 



In an infinite medium the integral of V • {pu}) vanishes as, by virtue of Gauss theorem, it can be reduced to a surface 
integral over a sphere with radius tending to infinity (40| . Therefore, if q(x, t) is known, the minimum of S is achieved 
at UJ = 0, and with p defined as explained above. As a result, Eq. (jAllI) yields the second equality in Eq. p3|) . 
Now we calculate the variation of S: 



SS 



dt / a-^ [ftu - D V (V • u)] • - DV {V ■ Su) - D' (V (V • u)) (V • ^u)] 



(A12) 



- ^[9tU-i5V(V-u)]'(V-(5u)| + y"dxJ-'((7(x,0))(5g(x,0) =0. 

Using Eq. (|A10p with a; = 0, we can rewrite the double integral in Eq. (|A12p as dt J d:ii}V with 

W = -Vp- [dtSu-DV{V ■ Su) - D' (V (V • u)) (V • 5u)] - — {Vpf (V • Su) . 

The first term, — Vp • dt6u, can be integrated over time by parts. Then, grouping the terms with D and D' together 
and employing Green's first identity and the boundary condition u(|x| — > oo, i) = 0, we obtain 



5S 



dt / dx 



dtp + DV^P+ -a' {Vpf (V-(5u) 



dx5u(x,0)- Vp(x,0) - VJ''(g(x,0)) 
J dx5u{x,T) ■ Vp(x,T) = 0. 



(A13) 



Demanding the bulk term to vanish we get Eq. the second of the MET equations. Each of the two boundary 
terms in Eq. (jAlBI) must vanish independently. The first term yields the boundary condition at < = 0: 

p(x, 0) - T'[q{x, 0)]=A = const. 

In view of the boundary conditions (jA9p and (jA4p and relation J-'{n) — 0, we obtain A = which yields Eq. ([TT|) . 

Now let us turn to the second boundary term in Eq. (jAlSP and assume that the domain is simply connected. 
As it turns out, the momentum field p{x,T) experiences a jump at the boundary 9f2 of the domain f2. To have a 
minimum action the third term in Eq. (|A13I) must vanish: 



J dx Su{x, T) ■ Vp(x, T) = . 



(A14) 
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This equality should be valid for an arbitrary (5u(x, T) compatible with the boundary conditions. In particular, it 
must vanish for any sufficiently smooth (5u(x, T) that vanishes at x S 51 U dfl. This means that Vp must vanish at 
X e ri, where = M"^ — is the exterior of the domain $7. This, along with the condition (|A9p . yields 

p(xen,r) = o. (A15) 

Hence Eq. (|A14[) can be rewritten as 

J dx 5u(x, T) • Vp(x, T) = 

n 

for any ft which contains fl and its boundary dfl. Integrating this equation by parts and taking into account Eq. (jA15p . 
we obtain: 

J dxp{x,T)S/- 5u(x, T) + J dS Sun (x, T) p„ (x, T) = , 

SI xeasi 

where Sun = 5u ■ n, n is the external unit normal to dfl, and Pini^,T)\_^^g^ = limx^-ao (xgO)P(x, T) is the value of 
p(x, T) along the boundary d^l inside fl. As Sq{x, T) = V • (5u(x, T) — Q inside il, the first term in Eq. (jA14l) vanishes. 
Now we need to deal with the second term: 

J rf5(5u„(x,T)p„(x,r) = 0. (A16) 

Let us apply Gauss theorem to the following auxiliary integral: 

J dSun{x,T)= J dxV-u(x,r). (A17) 

xedn xef2 

As q{x e ri, r) = k(x), the value of integral in the right side of Eq. (IA17P is fixed, see Eq. (jA3[) . Hence, the variation 
i5u must obey J^^g^dS Sun{x,T) = 0. Comparing this result with Eq. (|A16|) we see that pi„ must be constant along 
dfl. 

Overall, the complete set of boundary conditions at i = and t = T for the annealed setting reads: 

p(x,0)= J-'[9(x,0)]; 



g(xe r!,T) = k(x); p(xer!,T)^0; j^^^, = const, (A18) 

where the constant is unknown a priori and is a part of the solution. With these boundary conditions, the solution 
turns out to be quite simple, see Sec. IIII Al because of the presence of a local integral of motion p = J-'{q)- In 
particular, one can immediately see from this integral of motion that q(ic g f2,T) = n. Furthermore, in the particular 
case of void, g(x g 0,T) = 0, one obtains p(x e r2,T) = — oo for all lattice gases with J^'(O) — —oo. These include 
the RWs, the SSEP, the KMP, and the ZRP with a(0) = 0. 



2. Quenched setting 



In the quenched case (/(x, 0) = n is specified, and the second term in Eq. (|A7p vanishes. To proceed, we can simply 
put J" = in Eq. (|A13I) leading to 



6S 









dx 






(ix(5u(x 


0). 



pt + D\/^p+ -a' {\/pf (V-(5u) 



+ J dx Su{x, 0) • Vp(x, 0)- J dx 5u{x, T) ■ \/p{x, T) = 0. 



(A19) 



We see that the bulk term in (IA19P and the second boundary term in time again yield Eq. ^ and the boundary 
conditions (|A18p . respectively. Integrating the first boundary term by parts and assuming that p(x, 0) is sufficiently 
smooth, we see that the first boundary term vanishes because (5g(x, 0) = V • (5u(x, 0) = 0. 
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As one can notice, q{x e il, T) and p(x G Jl, T) are unknown a priori. For a class of lattice gas models (which include 
the RWs, the SSEP and the KMP, among others), one can show that, as in the annealed setting, p(x e ri,^) = — oo 
in the particular case of void at t = T. This a priori knowledge is convenient for solving the MFT equations, see 
Sees. [WlandlVl 

One can also notice that, by fixing (/(x, 0) = n, we narrow the class of density variations in the problem of minimizing 
S, in comparison with the annealed case. As a result, S in the quenched case is always greater than, or equal to, S 
in the annealed case. Therefore, T'qucnchcd < "Panncaicd- 



Appendix B: Microscopic theory for non-interacting random walkers 



In this Appendix we outline a microscopic theory of void formation for the RWs, both in the quenched and annealed 
settings. This approach yields exact results for the void formation probability and for the expected final density profile. 
In the long time limit, exact results coincide with predictions of the MFT formalism. For simplicity, we present the 
exact results only in the one-dimensional case. 



1. Quenched setting 

In the quenched setting, the initial condition is deterministic. We assume that, at f = 0, RWs occupy an infinite 
one-dimensional lattice, so that there are rim particles at site to. We demand that the [— L, L] interval of the lattice 
be empty at time T. Consider a single RW starting at t = at site to. This RW will be at site j with probability 
e~'^'^ Ij-m{2T), where Ik is the modified Bessel function (see, e.g., [25]). The probability for this RW to be outside 
the [—L, L] interval at time T is equal to the sum of arrival probabilities of this RW to all sites j such that |j| > L: 

e"'^ E Ij-rni^T) (Bl) 

\j\>L 

As all n„i RWs on the same site to are independent of each other, the probability that all of them will be found 
outside the [—L,L] interval at time T is 



-2T 



\3\>L 



(B2) 



Now, the RWs at other sites are also independent, so the probability P{L, T) of the void formation at time T is equal 
to an infinite product of probabilities (|B2I) : 



'^quenched (-^; T) 



n 



-oo<m<oo 



\j\>L 



(B3) 



This is an exact result, valid for any deterministic initial configuration. In particular, for a deterministically homoge- 
neous gas n„j = n = const, we obtain 



'^quenched (-^1 



n 



- oo < < oo 



-2T 



E I^-rnim 
\j\>L 



(B4) 



The most interesting asymptotic behavior arises when L ^ oo and T — cx), while the ratio £ = L/\/AT remains finite. 
In this scaling region we can employ the asymptotic relation 



e-^^/„(2T) 



1 



X 



which holds when T — >■ cx) and to — >■ cx3, while X is finite. This yields 



4T 



_2T V- , ^ erfc(^-X) + erfc(£ + X) 



(B5) 



(B6) 
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Equation (jB6|) allows one to recast Eq. (IB4[) into the scaling form (|40|) and pT|) . 

We can also use the exact microscopic approach to determine the average gas density outside the void a.t t = T. 
The probability that a RW starting at site m is at site x at time T is e~^^/2;-m(2T). This probability is obtained 
by sampling over all evolution histories, while we must consider a subset of histories when the RW ends up outside 



the [— i, L] interval. The corresponding conditional probability is e 
Therefore, the average density outside the void at t = T is given by 



2T 



I. 



-2T 



3\>L ^3- 



{2T) , see dHU). 



qi\x\>L,T) 



^ e"^^^'4_„(2r) 



(B7) 



In the limit of a; — >■ oo, T ^ oo and finite X = a;/V4T, we can use the asymptotic (jB5[) . As a result, the average gas 
density dX t = T acquires the scaling form (|49|) . 



2. Annealed setting 



In the annealed setting the initial state of the gas exhibits equilibrium fluctuations. In this case, the evaluation of 
the probability of the void formation at time T should be performed in two steps: 

1 . Evaluate the probability of the void formation at time T when starting from an arbitrary but fixed microscopic 
configuration (that is, for a quenched setting). 

2. Average this probability over the microscopic equilibrium density distribution. 



We have already performed step 1 of the calculation, and arrived at Eq. (|B3p . It is convenient to temporarily deal 
with a finite lattice which includes 27V + 1 sites: — iV, —N + 1, . . . ,m, . . . , N — 1, N, where N ^ L. The probability 



of void formation at time T when starting from the microscopic configuration (n_Ar, n^Ar+i, , 



at time i = is 



'^quenched (-^, 



n 



'N<m<N 



-2T 



\j\>L 



,nN-i,nN) 



(B8) 



We can now proceed to step 2. The (Poisson) microscopic equilibrium density distribution with a fixed total number 
of particles K can be represented as 



iV+l, 



. , Tim, . . • , nN-l,nN) 



where S is the Kronecker delta imposing the constraint of constant total number of particles: ^ 



is left is to average the quenched probability (|B8|) over the distribution (|B9 

"Paimoalcd = ^ '^quonchod(i, T) "^{ri-N, ■ ■ ■ , 



\m\<N "•™ 



(B9) 
K. What 

(BIO) 



The resulting expression is exact, but a more useful information can be extracted if we take the limits of T — !■ oo and 
X — >■ oo, keeping £ = L/ \/4r finite. Using Eq. (|B6|) . we obtain 



111 "^annealed — ^ ^ 

In 

\7n\<N 



erfc(£ -~X)+ erfc(£ + X) 



(Bll) 



instead of Um, we obtain 



where X = m/v 4T. Using a continuous spatial coordinate x instead of m and a coarse-grained gas density q{x, t = 0) 

POO 

InPannoaicd — / dx q{x , 0) In x{x , T) , (B12) 
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where xi^^T) = i[erfc(£ — X) +erfc(^ + X)] with X = x/\/4T, and we have returned to an infinite system. Similarly, 
we can go over to the continuous limit in the exact microscopic density distribution (|B9p . Employing Stirling's 
formula, ln(nm!) ~ Um In rim — and replacing sums over m by integrals over x, we obtain, up to a pre-exponent. 



This is simply the Boltzmann-Gibbs cost of the macroscopic density profile q{x, 0) in the RWs model, see Eq. 
and Table 1. Now Eq. (jBlOp becomes a path integral: 



which can be evaluated by the saddle-point method. The saddle point, in the functional space of q{x,0), is found at 
q{x,0) = nx{x,T) which coincides with Eq. (|25|) for the most likely initial density profile as predicted from the MET 
formalism. With this optimal initial density profile, Eq. (jB14|) yields the equilibrium result InT^anneaicd — —2nL. 
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